Missing derivative discontinuity of the exchange-correlation energy for attractive 

interactions: the charge Kondo effect 



(N 

o 

< 

in 



-a 

o 
o 



> 

O 

cn 

od 
O 
(N 



X 



E. Perfetto 1 and G. Stefanucci 1 ' 2 ' 3 

1 Dipartimento di Fisica, Universita di Roma Tor Vergata, 
Via de.Ua Ricerca Scientifica 1, 1-00133 Rome, Italy 
2 INFN, Laboratori Nazionali di Frascati, Via E. Fermi 40, 00044 Frascati, 
1 European Theoretical Spectroscopy Facility (ETSF) 



Italy 



We show that the energy functional of ensemble Density Functional Theory (DFT) [Perdew et al., 
Phys. Rev. Lett. 49, 1691 (1982)] in systems with attractive interactions is a convex function of 
the fractional particle number TV and is given by a series of straight lines joining a subset of ground- 
state energies. As a consequence the exchange-correlation (XC) potential is not discontinuous for 
all N. We highlight the importance of this exact result in the ensemble-DFT description of the 
negative-C/ Anderson model. In the atomic limit the discontinuity of the XC potential is missing 
for odd N while for finite hybridizations the discontinuity at even N is broadened. We demonstrate 
that the inclusion of these properties in any approximate XC potential is crucial to reproduce the 
characteristic signatures of the charge-Kondo effect in the conductance and charge susceptibility. 

PACS numbers: 71.15.Mb, 05.60.Gg, 31.15.E-, 72.10.Fk 



Density functional theory 1 ' 2 (DFT) provides a rigor- 
ous and computationally viable tool to calculate the elec- 
tronic properties of many-particle interacting systems. In 
spite of the great success in a wide range of applications, 
its practical use is still problematic in systems with fluc- 
tuating number of particles. Popular approximations like 
LDA or GGA are inadequate to predict, e.g., the band 
gap of solids, 3 ' 4 the correct dissociation of heteroatomic 
molecules 5-7 or the electrical conductivity of nanoscale 
junctions. 8 ' 9 A conceptual advance to deal with these 
cases is the ensemble-DFT put forward by Perdew et al.. 5 
These authors extended the original DFT formulation 1,2 
to a fractional number N of electrons, and pointed out 
the non-differentiability of the energy functional E[n] of 
the density n at integers N = J n. Typically the discon- 
tinuity in dE/dN is the difference between the ionization 
energy and the electron affinity since for any N between 
two consecutive integers M and M + 1 one has 



E(N) = (M + 1 - N)£ M + (N- M)£ m +i 



(1) 



i.e., E(N) is a series of straight lines joining consecutive 
ground-state energies £m of the isolated system with M 
particles. Fig. 1 (top panel) illustrates the typical out- 
come of a ground state calculation of E(N). It is worth 
recalling that the crucial hypothesis for the validity of 
Eq. (1) is the convexity inequality 
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£m+i + £m-i — 2£m > 0. 



(2) 



Indeed, in this case one can show that the density matrix 
which minimizes the total energy is a linear combination 
of projection operators over the ground states with M 
and M + 1 particles. As discussed in Ref. 5 the hy- 
pothesis Am > is certainly reasonable in systems with 
repulsive interactions. In contrast the convexity inequal- 
ity can be violated in the attractive case. For instance 
in the attractive Hubbard model Am is positive for even 
M and negative otherwise. 10,11 What consequences docs 
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FIG. 1: Illustrative examples of the total energy E(N) (solid 
line) as a function of the fractional number iV when the con- 
vexity inequality (2) is (top panel) and is not (bottom panel) 
fulfilled. In the bottom panel the dashed line joins consecutive 
ground-state energies and differs from E(N). 



the break-down of Eq. (2) have in ensemble-DFT? What 
are the physical implications? 

In this paper we generalize Eq. (1) to arbitrary Am- 
In particular we prove that E(N) is a convex function 
given by a series of straight lines joining a subset of 
ground-state energies, as schematically illustrated in Fig. 
1 (bottom panel). We further study the implications 
of the missing derivative discontinuity in the negative- U 
Anderson model. In this prototype system the attrac- 
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tive interaction is at the origin of the so called chargc- 
Kondo effect: 12,13 at very low temperature the fluctua- 
tions between the empty state and the doubly-occupied 
state of the impurity level produce a strong enhancement 
of the charge susceptibility, accompanied by a drastic 
shrinkage of the conductance peak. We show that the 
key features of the charge-Kondo effect can be captured 
within cnscmblc-DFT. We extend a recently proposed 
functional 14 devised for the spin-Kondo effect to attrac- 
tive interactions and account for the broadening of the 
discontinuity 15 due to the finite hybridization of the im- 
purity level. The transition from the spin-Kondo effect 
to the charge-Kondo effect is caused by the shift of the 
discontinuity of the exchange-correlation (XC) potential 
from N = 1 at U > 16,17 to N = and N = 2 at U < 0. 

Theorem : Given the ground-state energies {£/} of the 
isolated system with I particles, if 



£i — Em £p — £m oj — £m 
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P-M 



J-M 



(3) 



for every I < M and every J > M, then in the range 
M < N < P it holds 



, . P-N N-M 



(4) 



Graphically this means that for N G [M, P] the energy 
E(N) lies on the straight line connecting Em to Ep if 
and only if the slope £ pZ_\f is larger than all the slopes 
of the lines connecting Em to £i<m an d smaller than all 
the slopes of the lines connecting Em to £j>a/, see Fig. 
1 bottom panel. 

Proof: We have to show that in the range M < N < P 

^l^lEl cannot be 
(4) for any {uj L } 



the variational energy E var (N) 
smaller than the energy E(N) in Eq 
constrained to satisfy 



and u>l > for all L. Using Eq. (3) one has 
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Y J U L (L-M) = E(N), (6) 



which proves the theorem. 

Thus E(N) in Eq. (4) is a convex function of N and 
reduces to Eq. (1) provided that the convexity inequality 
Am > is satisfied for all M since in this case P = M+l. 
The physical content of the theorem is clear. If the sys- 
tem is open to a charge reservoir the density matrix at 
zero temperature is a mixture of ground states with M 
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FIG. 2: E(N) of Eq. (4) (solid line) for the 4-sites attractive 
Hubbard ring with U = —5 and on-site energy Ed = 3.5. 
Energies are in units of the hopping integral. The dashed line 
joins consecutive ground-state energies. 



and P particles. For example in Fig. 2 we show that in 
the attractive Hubbard ring P = M + 2, see also Refs. 
10,11. The value P = M + 2 is peculiar of attractive 
systems where the electron pairing causes Am ^ for 
even/odd M. This property is consistent with the exper- 
imental observation of the Coulomb blockade of Cooper 
pairs 18,19 (Cooper staircase) in superconducting single- 
electron transistors, where a superconductive island is 
connected to metallic leads. In this situation the appli- 
cation of a gate voltage v g to the attractive region causes 
a jump of 2 in the number of particles at the special 
values v g = ^£z^Jf±a .ao 

In ensemble-DFT the discontinuity dE/dN is the sum 
the Kohn-Sham (KS) discontinuity, which is zero for odd 
N, and the XC discontinuity A xc (iV). Since dE/dN = 
for odd N we conclude that 



A KC (N) = 



for odd N. 



(7) 



In the following, we consider a negative-t/ Anderson 
model as an example in which XC discontinuity is miss- 
ing. The Hamiltonian describes a set of noninteracting 
electrons coupled to a site at which Hubbard type interac- 
tion occurs. 12 ' 13 This is an effective model for conduction 
electrons coupled to an interacting impurity with vibra- 
tional modes. For strong electron-phonon coupling the 
polaronic shift can overcome the Coulomb charging en- 
ergy and the effective electron-electron interaction turns 
out to be attractive. The Hamiltonian reads, in standard 
notation, 



Q,cT i—1 <J 

+ Un dt n di + t'Y^[c\aa d v + h -C-}, (8) 



where t is the nearest-neighbour hopping in the leads, 
t' is the lead-impurity hopping, U < is the attractive 
interaction, and v g is the gate voltage coupled to the 
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FIG. 3: Potential « Hxc at very low temperature, /3 = 400/|[/|, 
and comparison between t!g xc at J3 = 1/W-y ~ 25/|£/| and 
uh xc at zero temperature with 7 = 0.125|C/|. The Lorentzian 
broadening is much larger than the thermal broadening for 
7 < \U\. 



impurity density n da — dj r d a . Below we focus on the 
half-filled system and hence take the chemical potential 
/i = 0. At very low temperature and gate voltage around 
v g = — U/2 = \U\/2 this model exhibits the so called 
charge-Kondo effect. 12,13 This effect consists in the for- 
mation of a "local pair" at the impurity and is due to 
strong charge fluctuations between the nearly degener- 
ate states |0) and | tl) of the empty and doubly occu- 
pied d-level. As predicted by Taraphder and Coleman 13 
the local pair is "screened" by the surrounding conduc- 
tion electrons and forms an "isospin singlet" . With in- 
creasing \U\ the main features of the charge-Kondo ef- 
fect arc i) the shrinkage of the conductance resonance 
at v f , = v„ and ii) the large growth of the charge sus- where 



ceptibility \d = —dn d /dvg. 13 ' 21 These results can be 
qualitatively understood by mapping the Hamiltonian 
of Eq. (8) into the positive-?/ Anderson model. Under 
a particle-hole transformation in the spin down sector, 
d\, — > d\ and Ciai — > (— l) l cj a , , the original Hamiltonian 
is transformed into the positive-?/ Anderson Hamiltonian 
with fixed gate voltage —\U\/2 and effective magnetic 
field B c g = —\U\/2 + v g coupled to {n#[ ~ n di _) . 22 ~ 24 
Since the magnetic field suppresses very efficiently the 
Kondo correlations 21-23, 25-28 the spin-Kondo effect in the 
transformed Hamiltonian occurs only in the proximity of 
Vg = Vg. Consequently the conductance drops rapidly to 
zero as v g deviates from v g . At resonance the spin fluc- 
tuations in the transformed Hamiltonian correspond to 
"isospin", i.e., charge, fluctuations in the original Hamil- 
tonian, thus leading to the formation of an isospin sin- 
glet (local pair). This phenomenology explains the large 
growth of the charge susceptibility Xd as \U\ increases 
(this growth is not observed for positive U). 

Let us show how these features can be captured in 
ensemble-DFT. In a recent Letter 14 an approximate 
Hartree-XC potential for the positive-?/ Anderson model 
was proposed. The exact energy functional of the isolated 
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FIG. 4: Conductance G as a function of v g for two values 
of U = —2, —6 (in units of 7) and different approximation 
schemes. The data is compared to the variational results of 
Ref. 21, which agree closely with the NRG data of Ref. 26, 
and can be, therefore, considered as exact. 



impurity reads 



'Axc( n d) = j+g{n d - 1), 



. , U la; 
9{x) = — + -x log - 



+ \R 



1 



(9) 



(10) 



(3 being the inverse temperature. For U > and in 
the limit j3 — > 00 the potential v^ KC (nd) — > U8(nd — 1) 
which has a discontinuity U at = 1. In the Wide 
Band Limit Approximation (WBLA) t, t' ~ > 00 with con- 
stant 2t' 2 /t = 7 -c U (weak tunneling rate) one can 
approximate the Hartree-XC potential on the impurity 
vhxc ~ «h xc and set it to zero in the leads. 14 The discon- 
tinuity forces the occupation to be unity for gate voltages 
< v g < U. 17,29 Thus the KS potential is pinned at the 
Fermi energy and the KS conductance exhibits a Kondo 
plateau as a function of w ff . 14 ' 30,31 

The physical argument leading to Eq. (9) is inde- 
pendent of the sign of U and we may argue that the 
functional i>h xc should predict, at least qualitatively, the 
correct conductance also for negative U. In the anal- 
ysis below we consider the zero-temperature case. For 
U < the potential Wy xc is not discontinuous at = 1 
but instead develops two discontinuities (of size |?7|/2) at 
n d = and n d = 2, see Fig. 3. 32 Within the WBLA we 
determine the occupancy on the impurity by solving the 



4 



14 
















u=o 




12 








U=-3 












— U=-6 




i n 








— U=-9 











o 

Xd 








6 








4 








2 


^ J- 










-0.4 -0.2 0.0 0.2 0.4 

(v s +U/2)/y 



FIG. 5: Charge susceptibility at the impurity \d = —dnd/dvg 
as a function of v g with potential Uh xc for different values of 
U (in units of 7). 



In Fig. 3 we show the comparison between v^ xc at finite 
temperature and ^jj xc at zero temperature. Choosing 
(3 = I/W7 = tt\U\/(8j 2 ) we see that the thermal broad- 
ening is much smaller than the Lorentzian broadening for 
7 -C \U\. Figure 4 clearly illustrates the crucial role of 
the broadening of the discontinuity in the shrinkage of 
the conductance resonance as \U\ increases. The figure 
displays also the conductance in the Hartrcc-Fock (HF) 
approximation, i.e., with potential ^HF^d) = Und/2. 
Even though this potential reproduces the shrinkage up 
to U ~ —2, it becomes unreliable already for U ~ —3. 
At this critical value the self-consistent equation for the 
density develops multiple solutions, three in our case, as 
shown in the bottom panel of Fig. 4. This multistability 
scenario should be contrasted with the positive-C/ Ander- 
son model where multiple solutions within the Hartrec- 
Fock approximation are found only out-of-equilibrium. 33 



self-consistent equation rid = 1 — f tan 1 " a+, ' I ! y xc 

with «Hxc = ^Hxc O nce rid is known we calculate the KS 
conductance from 



Go [v g + VHxc(«-d)] 2 + 7 2 ' 

where Go = 1/V is the quantum of conductance. The 
exact KS conductance equals the exact conductance due 
to the Friedel sum rule and the WBLA. 9 It can easily be 
seen that for wh X c = u h xc the conductance is correctly 
peaked at v g = v g but its width is weakly dependent 
on U. Indeed V^Jjid) — U/2 everywhere except that 
at the occupations rid = 0, 2, see Fig. 3. Therefore the 
conductance as a function of v g has a constant width 7 
since rid is never exactly or 2. This is illustrated in 
Fig. 4 where the conductance calculated using -y^ xc is 
compared with the variational results of Ref. 21. 

The potential w^ xc can be substantially improved by 
following the observation of Ref. 15. At temperatures 
T = 1/(3 below the Kondo temperature the broadening 
of the discontinuity in i^ xc is proportional to pUe~^ u ^ 2 
and approaches zero in the limit T — > 0. However, the ex- 
act Hartrec-XC potential should have an intrinsic broad- 
ening W ~ 87/71- 1 [/ 1 due to the finite hybridization of the 
impurity. Therefore we here propose a Hartree-XC po- 
tential which is the convolution of v^ xc with a Lorenzian 
of width W. The resulting potential for negative U and 
zero temperature reads 

1 , \ U U . nd — 2 nd . , 

«HxcM = — + — (arctan — + arctan — ) . (12) 



Finally we used the Hartrec-XC potential v Kxc to cal- 
culate the charge susceptibility \d = —dnd/dv g . In Fig. 
5 we show Xd as a function of v g for several values of U. 
Also in this case our approximation correctly captures 
the growth of \d at v g = v g , another typical signature 
of the charge-Kondo effect. We further observe that the 
hight of the peak in Xd saturates to values around 2 if 
we use the Hartree-XC potential t>^ xc (not shown). Thus 
the broadening of the discontinutity is crucial in this case 
as well. 



In conclusion wc generalized the variational energy 
functional of ensemblc-DFT to cases where the convexity 
inequality is not fulfilled. The energy E(N) is a convex 
function of the fractional particle number N, and it is 
given by the lowest series of straight lines joining a sub- 
set of ground-state energies. We discussed the relevance 
of this property in the description of correlated systems 
with attractive interactions. As for odd N the energy 
E(N) has no cusp and the KS discontinuity is zero, the 
XC discontinuity is zero in these cases. We showed that 
the missing XC discontinuity and the broadening induced 
by the finite hybridization with the leads are essential 
features of any approximate functional to describe the 
charge-Kondo effect in the negative-?/ Anderson model 
within ensemble-DFT. The functional proposed in this 
work yields results in fairly good agreement with the 
available numerical data. In particular the shrinkage of 
the conductance peak as well as the growth of the charge 
susceptibility with increasing \U\ are correctly captured. 
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